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I.  INTRODUCTION 


Occasionally,  the  need  arise i  for  the  analysis  of  moderately  thick 
shells  with  rapidly  varying  wall  ihickness.  Local  stress  concentration 
becomes  the  priti.ary  concern  in  certain  oases  rather  than  the  averaged 
behavior  approx Imated  by  conventional  thin  shell  analysis.  Clearly  the 
determination  of  such  information  requires  a  careful  consideration  of  the 
equations  of  tliree-dimensional  elasticity  theory. 

In  this  repcrt,  we  initially  consider  the  class  of  moderately  thiok 
shells  of  revolution  under  axially  symmetric  loadings,  i.e.,  the  deformation 
will  be  axially  symmetric.  The  local  effects  of  geometrical  irregularities 
in  the  wall  thickness  (e.g. ,  circumferential  grooves  or  nodes)  will  be  of 
principal  concern.  Although  the  corresponding  problem  class  in  three- 
dimensional  elasticity  is  well-known,  for  example  Reference  1,  exact 
solutions  in  general  are  extremely  difficult  to  obtain.  The  plan  of  analysis 
therefore  consists  in  utilizing  the  relative  thinness  of  the  shell  to 
systematically  approximate  the  three-dimensional  equations. 

A  plan  of  successive  approximations,  conceptually  due  to  0.  Gdhner,  * 
is  adopted.  The  choice  of  the  order  of  approximation  is  obviously  sensitive 
to  the  loadings  considered.  A  particular  choice  is  considered  in  this 
report  leading  to  a  structure  permitting  the  use  of  plane  biharmonic 
functions. 

It  is  shown  that  this  plan  is  suitable  for  handling  a  restricted  but 
useful  class  of  loadings  of  cylindrical  shells  with  rapidly  varying  thick¬ 
ness.  In  particular,  the  zero'th  order  solution  supports  Neuber’s  assertion* 
that  when  the  shell  is  very  thin  for  a  class  of  applied  loadings  the  problem 
can  be  considered  as  one  of  plane  strain.  The  first  two  orders  of  solution 
oan  be  explicitly  formulated  in  the  framework  of  analytic  functions  of  a 
complex  variable. 


II.  FORMULATION  IN  TERMS  OF  SUCCESSIVE  APPROXIMATIONS 

For  the  problem  class  of  axially  symmetric  deformation,  it  will  be 
convenient  initially  to  review  the  formulation  in  terms  of  cylindrical 
coordinates  (r,  6,  z)  where  z  is  the  axis  of  revolution.  Assuming  no  body 
forces,  it  is  well-known,  e.g.,  Reference  1,  that  the  equations  of  equilibrium 

(1) 

Tr,,r  ♦  *  Tr./r  "  ° 

where  an  independent  variable  appearing  as  a  subscript  and  following  a 
comma  denotes  partial  differentiation  with  respect  to  that  variable.  The 
compatibility  equations  for  this  case  are 


-3- 


(2) 


where 


V*<rf  -  2(<rf-  <r$)/ r*  ♦  n,  r r/(l  ♦  v)  •  0, 

♦  2(<rf  -  ♦  n,f/r(l  ♦  v)  ■  0, 

V*o-f  ♦  ♦  v)  •  0, 

V*Tr.  "  Tr*/r*  ♦  ^.,g/(l  ♦  *)  ■  0; 

V*  ■  aVar*  ♦  (l/r)3/3r  ♦ 


n  ■  (a,  ♦  <r,  ♦  a§ )/ ( 1  ♦  v), 
and  v  denotes  Poisson's  ratio. 

The  systems  (l)  and  (2)  can  be  replaced  by  a  formulation  in  terms  of 
a  single  stress  function,  <f>,  as  proposed  by  Love,4  by  defining 


<tt  -  [i JQl<t>  -  <£„]., 

<r$  -  [v7*d>  -  4>t/r]t, 

(3) 

»■,  *  t(2  -  v)Vl4>  - 

Tr«  * 

where 


V*V*0  ■  0. 


(in  (3),  subscript  notation  denotes  partial  differentation  of  the  stress 
function  <t>(r,  z ) . ) 

Let  us  now  introduce  the  new  coordinates  £,  rj  defined  by 

£  -  r  -  R,  tj  -  z,  (4) 

where  R  is  an  appropriately  chosen  fixed  radius  of  the  shell.  For  the 
present  problem  class,  i.e.,  cylinders  of  variable  thickness,  we  consider 
R  as  independent  of  r  and  z.  The  previous  Bystems  can  now  be  considered 
in  terms  of  the  £,  tj  coordinates.  For  example,  the  equations  of  equili¬ 
brium  and  compatibility  have  the  forms 

°>,f  *  T(v,v  *  {<Tf  *  "  °’  (1'} 
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(2*) 


where 


t*e  ♦  *  R)  -  2{oy  -  <t9)/[€  ♦  R)*  ♦  ♦  v)  -0, 

a  ■  a*/d£*  ♦  3*/V. 


Gfihner's  plan  of  successive  approximations  depends  on  the  validity  of 
the  expressions 


l/(f  ♦  R)  ■  (1/R)  f  <-  f/R)\ 
l/(f  ♦  R)1  -  (l/R)1  Z  (nd){-  £/R)>. 


(5) 


Ve  shall  henceforth  confine  our  attention  to  regions  where  the  expansions 
(&)  can  be  considered  valid.  Thus,  the  expansions  (5)  will  be  considered 
in  (l*)  and  (2*)  and  a  plan  of  successive  approximations  to  the  system 
will  be  defined.  In  particular,  if  a  is  the  generic  symbol  for  the 
stresses,  it  is  assumed  that 


a 


2  o'*1 
*■0 


* 


where  the  superscript  denotes  the  order  of  the  approximation. 


(6) 


In  the  method  of  successive  approximations,  the  plan  of  approximation 
is  dearly  not  unique.  Furthermore,  a  rigorous  assessment  of  the  validity 
of  any  given  plan  is  generally  unfeasible.  On  the  other  hand,  several 
minimal  requirements  are  obvious.  The  plan  should  lead  to  a  systematic 
satisfaction  of  equilibrium  and  compatibility  at  least  in  an  asymptotic 
sense.  Furthermore,  each  order  of  solution  should  oorrespond  to  a  well- 
defined  mathematical  system  of  equations.  Finally,  a  reasonably  accurate 
approximation  of  the  prescribed  load  conditions  should  be  possible. 


Ve  propose  the  following  particular  plan  of  successive  approximation!: 
For  the  n'th  order  approximation  with  leading  components  cr'1’,  terms  con¬ 
taining  £*/R,+1  explicitly  in  (l* )  and  [•&')  will  be  neglected.  Thus,  the 
sero'th  order  approximation  corresponds  to  the  system 


(0) 

(6) 

.  («»  m 

♦  T .  " 

tfn>  v 

°’  re r,.  0, 

A**’ 

n(9t  ,, 

♦  n.M/(i 

♦  v)  "0,  Cxr9*]  -  o, 

.  ♦  v) 

to 

10)  ,, 

♦  ♦  v)  -  o. 

(7) 
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For  n  ■  1,  the  corresponding  system  is 


1 1 1  Ill  .  (01  i#» ,  . 

T<v.v *  (oV  '  ae  >/R  "  °* 

in  in  itt  . 

T(v.<  *  Ttr,  /R  -  0. 

etc. 


(8) 


The  plan  oan  be  conveniently  eunmariaed  in  terms  of  stress  functions. 
Guided  by  the  structure  of  the  stress  function  (3),  it  is  a  straightforward 
matter  to  determine  a  system  of  stress  functions  compatible  with  the  pre¬ 
ceding  plan  of  successive  approximation.  In  fact,  the  n'th  order  approxima¬ 
tion  can  be  replaced  by 

i*  i  re.  „in  .  ,  .  »  ij-ii , 

<T(  -  JvAF  -  +  (v/n)  2Y(  (-  f/R) 

e»  f  .„lil  ,  ,  »  '} 

ort  •  |vAF  -  (1  -  w)  2^  (-  f/R)  /RJ„* 

e.  r  ...Hi  „(•>  ,  ,  •  «j-n,  , *-j  .1 

%  ■  ul  -  vW  *  Fi(  ♦  (2  -  v^^Fj  (-  f/R)  /®J,. 

r(v  *  '  vAF  ♦  <1  -  *0  (-  f/R)  /RJ^» 


O) 


n  ■  0,1, £,.... 


where  F 11  (£,17)  "0  if  i<0  and 


in 


MF’"  ♦  (E/Rj^'AFi"1'1’!-  f/R)‘ 

i  * 


(10) 


a ,  «-> 


»/ »-i 


.(■-1-SI 


-(i/R-)-z’(iti)F;:  (-f/R)  ♦(i/sr  )  2  (i*8)(i*i)p;  (-e/*)  -o 

i«t  *•  4*1  * 


n  •  0, 1, 2, .... 

The  plan  of  approximation  is  not  yet  complete.  Although  a  procedure 
has  been  invented  for  the  systematic  consideration  of  equilibrium  and 
compatibility,  there  is  no  assurance  that  each  system  is  well-defined  and 
that  reasonable  approximation  of  the  applied  loads  can  be  realised.  In 
order  to  study  the  oomplete  structure  of  each  Bystem,  it  is  necessary  to 
oonsider  the  matter  of  boundary  conditions  to  be  imposed  on  the  individual 
ayatem8.  Therefore  the  applied  load  or  possible  systematic  approximations 
of  the  applied  load  must  be  investigated. 
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III.  THE  ZERO' TH  ORDER  THEORY 


For  n  ■  O,' (9)  and  (10)  yield 


where 


r  id 

(til 

LvAF  -F«Jv 

r  t«> 

(#11 

-  V\a( 

(1  -  v)AF 

Itl  Itl 

+  F« 

r  i») 

,  ml 

■ 

** 

i 

vAF 

(0) 

MF 

0. 

(U) 


(12) 


Although  this  system  is  apparently  consistent  with  the  well-known  problem 
of  plane  strain,  such  a  conclusion  cannot  be  immediately  arrived  at  with¬ 
out  further  consideration  of  the  applied  load  conditions. 

In  Figure  1,  the  force  resultants  T^  ds  and  ds  per  unit  of  circum¬ 
ference  acting  on  the  arc  ds  ia  shown  where  it  is  assumed  that  the  foroe 
is  exerted  by  the  material  to  the  left  on  the  material  to  the  right  pro¬ 
ceeding  from  A'  to  B'. 


From  equilibrium, 

“  %  sih/3  ♦  T(n  008/8  «  df/ds  +  TM  irj/ da, 

T'f  •  oy  cos/3  ♦  r(r)  aiiv3  ■  cr^  dfj/ds  -  dtf/ds. 

We  now  define  as  an  additional  part  of  our  plan  of  successive 
approximations, 


■  ( j  i 

T.  T  ,  T* 


I  T^\ 


J-0 


where 


I*  l 
v 

id 


(■I  <•>  , 

dsyds  ♦  r(l)  d7}/ds,  n  -  0,1,2,. 


v 

(it 


Thus, 


lS  -  Vf  d-rj/da  -  rf  d£/ds,  n  ■  0,1,2,  ...  . 

T^**  -  -  d  F^’/ds  -  (l  -  df/ds  -  AF^**  dr>/ds^, 

lot  f  lot  .  101*1  , 

-  -  d|F^  -  v/SF  J/ds. 


(15) 

(14) 


(15) 


(15) 


(17) 


(18) 


Consider  now  the  resultant  forces,  R^  and  R*  acting  on  the  surface 
generated  by  the  arc  A'B'.  Rtf  vanishes  by  symmetry,  whereas 


R,  *  2 v  f*[R  *  f)  ds. 


(19) 


If  the  arc  A'B'  is  a  closed  contour  with  respect  to  the  frj-plane,  then 
should  vanish  as  B*  coincides  with  A'  from  a  consideration  of  equilibrium. 
If  T^  is  approximated  by  T1*1  and  substituted  in  (19),  it  is  easily  shown 
that  R  doe 8  not  vanish  in^general  for  a  closed  circuit.  It  is  therefore 
clear  that  a  meaningful  definition  of  R  must  be  given  to  each  order  of 
theory  allowing  for  the  order  of  approximation  to  the  equilibrium  state. 

A  similar  argument  can  be  made  for  the  displacements  and  the  resultant 
moments. 


A  well-defined  zero’th  order  theory  will  now  be  considered.  We 
assume  the  stresses  of  the  zero'th  order  theory  are  defined  by  (ll)  and 
(12).  The  function  F1*1^,^)  and  its  fourth  order  partial  derivatives  are 
therefore  assumed  to  exist  at  all  interior  points  of  the  region  described 
in  the  fi7-,plane.  Consider 


(30) 
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where  T'*'  and  T^*’  are  defined  by  ( 17)  and  (18).  For  any  closed  contour 
A*B*  (bounding  a  simply  connected  region  of  the  material  in  the  £i>-plane), 
P** 1  and  P**’  vanish  due  to  the  continuity  of  F***  and  its  third  order 
partial  derivatives.  Thus,  ( 11 )  and  (12)  coupled  with  (17)  and  (18)  imply 
a  system  of  applied  loads  whose  resultant  force  vanishes  from  a  plane  point 
of  view.  A  similar  result  follows  when  the  resultant  moment  of  the  applied 
load  is  considered. 

It  has  thus  been  shown  that  for  the  plan  of  successive  approximation 
outlined  above,  the  zero'th  order  approximation  implies  a  system  of  applied 
loads  which  is  telf -equilibrating  in  the  plane  eenee.  It  would  be  in¬ 
consistent  therefore  to  attempt  to  use  this  particular  plan  of  approximation 
for  loading  systems  whioh  violate  this  condition  to  the  first  order,  e.g. , 
internal  pressure  acting  on  a  cylindrical  shell.  Therefore,  we  henceforth 
restrict  the  class  of  applied  loads  to  those  which  are  self-equilibrating 
in  the  plane  sense  to  the  first  order. 


Several  details  are  still  necessary  to  correlate  the  zero'th  order 
theory  with  the  original  three-dimensional  problem.  The  foroe  resultant 
defined  by  (l9)  will  now  be  considered  as 


R 


n 


£  R 

!-• 


* 


(22) 


Therefore,  a  compatible  definition  of  R^**  with  the  zero'th  order  theory  is 


R^*'  ■  ZnR  T^*’  ds.  (23) 

Finally,  we  assume  the  radial  and  axial  displacement  components  u  and 
w  can  be  expressed  as 


s  lit  s  <1> 

£  u  ,  w  ■  £  w 

J ••  J"* 


(•24) 


It  will  be  recalled  that  for  the  zero'th  order  theory, 


m  .  (*)  (•). 

O-.  m  v[<r.  ♦  <r  ). 


(20) 


Therefore,  from  Hooke's  law  and  the  well-known  strain-displaoemsnt  relations 
in  cylindrical  coordinates, 
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E  3u  /3£  •  (l  -  v  )er(  -  v(l  ♦  , 

_  .  <•>  ...  .  i.  <•>  .  .  <•» 

E  3w  /3t>  •  (1  -  v  )<r^  -  v(l  ♦  v)a#  , 

0(3u  /3t?  *  3w  /3f)  -  rfv  , 


where  B  is  Young's  modulus  and  G  -  E/?(l  ♦  v ).  Thus, 

ir,#» 
r«l' 


2G  3u  /3£  ■ 


2G  3w  /3tj  ■  (l  -  gvjAF^  ♦ 


Apart  from  rigid  body  motions, 


2G  u 


<•> 


?ev  ' 


«•»  ,  „  ,A  «*) 

2G  w  -  (l  -  2v)AF  ♦  ¥(e  . 


(26) 


(27) 


(22) 


IV.  THE  ORDER  OP  THEORY,  n  -  1 

For  orders  of  theory  n>l,  the  stress  functions  F(>>(£,i?)  are  no  longer 
plane  biharmonic  in  general.  On  the  other  hand,  the  application  of  plane 
biharmonic  functions  can  be  preserved  provided  a  particular  integral  of  (10) 
is  found  for  each  order  of  theory.  In  general,  the  determination  of  the 
particular  integrals  requires  a  knowledge  of  the  preceding  orders  of 
solution. 


Fortunately,  the  particular  integral  for  the  order  n  •  1  can  be  found 
direotly.  In  fact,  for  n-  l,  (-f/R)F’#,/2  is  an  integral  of  (10).  There¬ 
fore,  the  function  G(1>  (£,t?)  is  defined  as 


Gm-  £F,#,/2R. 


(29) 
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Then, 


Ul 


f  4 in  _._isi .  _  _is»,_l 

£m4G  —  £  •  i^AF  /2R  ♦  fF^  /2R  ♦  /Rj^i 


III  r 

r#  d> 


4/,m  „  <•»  .  „ 

vbG  -  vfAF  /2R 


111 


;*'4. 

[(1  -  v)AG,I’+  o“’  -  (1  -  v)£AFm/2R  -  ^J'/2r],, 

r.in  -in  _<*i  .  „ 

[G^  -  vAG  -  fF^  /2R  ♦  wfAF  /2Rjj, 


(11 

rev 


end 


M  G,n  -  0. 


The  corresponding  stress  resultants  are 


ill  (ll...  ill 

T„  «  -  ar_ 

7 ) 


v  di/da  ♦  Tfv  dTj/ds 
dG^* '/da  -  (l  -  vJ^AGy'df/ds  -  AG^'dTj/ds^ 

-  fr^"/2R  -<[f“'  -  vAF"']dV®ds, 

T^1’  ■  ^‘Wds  -  r^'d^/da 

-  -  d[0;;’ .  vAG,n]/ds .  fr;#’/ER 
♦  F^'dij/Rds  ♦  [f;;*  -  I/AFl"]df/2Rds. 


(30) 


(31) 


(32) 


(33) 


An  interesting  mathematical  dilemma  occurs  from  a  consideration  of  the 
determinancy  of  the  solution.  In  general,  a  logioal  definition  of  T*1*  and 
T^11  would  lead  to  a  plane  self-equilibrating  set  of  prescribed  loads.  For 
example,  if  T^**  ■  T.  and  T^**  ■  on  the  boundary,  the  natural  choice 
for  the  boundary  conditions  ror  the  order  of  theory  n  >  1  would  be 
T***  »  Ti1’  ■  0  which  corresponds,  of  course,  to  a  plane  self-equilibrating 
system  or  applied  loads.  On  the  other  hand,  it  is  a  simple  matter  to  show 
that  in  general  T'1’  and  Ti 1  *  as  defined  by  (32)  and  (33)  are  not  self- 
equilibrating  in  the  plane  sense.  Although  the  portions  of  (32)  and  (33) 
involving  Q,*,(£»i|)  correspond  to  a  self-equilbrating  load  system  in  the 
plane  sense  by  arguments  similar  to  those  of  Section  3,  it  oan  easily  be 
Shown  that  the  remaining  portions  involving  the  sero'th  order  solution  in 
general  violate  the  conditions  of  plane  equilibrium.  Since  the  flexibility 
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of  the  solution  rests  in  the  determination  of  Gm(f,i?),  clearly  from  a 
precise  standpoint,  G*1^,^)  is  indeterminate. 

The  difficulty  described  above  is  common  to  most  procedures  which 
depend  on  systematic  approximation  of  a  system.  The  dilemma  can  be  by¬ 
passed  by  introducing  an  additional  approximation.  If  G1 1  *  (£,i?)  is  de¬ 
termined  on  the  basis  of  ensuring  prescribed  loads  over  a  portion  of  the 
boundary  only,  in  many  practical  problems  it  oan  be  argued  that  sufficient 
aoouraoy  of  solution  in  these  regions  is  assured  if  the  remaining  portions 
of  the  boundary  are  handled  in  only  an  average  sense. 

To  illustrate  the  argument  above,  wo  consider  edge  loading  of  the  shell 
with  load-free  conditions  on  the  lateral  surfaces.  It  is  a  straightforward 
matter  to  show  that  T'1'  ♦  (^/R)T ‘ • *  and  T^1’  ♦  (f/R)T^*>  correspond  to 
plane  self-equilibrating  load  systems.  On  the  lateral  surfaces  of  the  shell 
T***  and  T^*’  vanish,  and  the  conditions  on  G,n  (f,rj)  on  the  lateral 
surfaces  can  be  found  from  (32)  and  (33).  On  the  ends,  Gm  (£,y)  oan  be 
adjusted  to  satisfy  the  edge  loading  in  an  average  sense.  By  appealing 
to  Saint  Venant's  principle,,  this  approximation  should  still  yield  a 
reliable  estimate  of  the  lbeal:  behavior  ini  the  regions  of  rapidly  varying 
thickness  provided'  these  regions  are;  sufficiently  removed  from  the  edges. 

(|  |)l 

Hie  arlBl1  force '  resultant'  R  consistent  with  (22)  for  this  order 

theory/  ib  defined1  as' 

R*1’  -  2rrR  £  [t‘“  ♦  (f/R)T^#,]ds.  (34) 

This  definition  of  in  addition  to  being  consistent  in  form  with  our 

plan  of  successive  approximation,  enjoys  the  property  of  defining  the  same 
total  resultant  axial  force  independent  of  the  cross  section  chosen.  This 
last  property  requires  that  (34)  be  independent  of  path  when  considered  a 
line  integral.  It  is  a  straightforward'  matter  to  show-  that  (34)  meets 
these  requirements. 

The  corresponding  displacement's  u(1>>>  and'  w'11’  oan  be  carried  out  in  a 
direot  manner;  however,  the  details  will  not  be  included  here. 


V.  REMARKS  ON  ORDERS  OF  THEORY,  n>l 

For  n>l,  the  approximations  can  be  analysed  in  a  manner  similar  to 
Sections  III  and  IV.  Again  the  problem  for  each  order  can  be  resolved  into 
the  determination  of  a  biharmonic  stress  function  along  with  a  particular 
solution  of  (10)  in  terms  of  the  preceding  orders  of  solution.  The  deter¬ 
mination  of  the  particular  solutions  becomes  increasingly  difficult  with 
increasing  n.  For  example,  for  n  ■  2  the  particular  solution  requires  a 
solution, of  the  Poisson  type  differential  equation  with  a  right hand  side 
proportional  to  F,,( 
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The  indeterminacy  of  eaoh  of  the  higher  orders  of  solution  rsseMbles 
the  ease  n  ■  1.  Thus,  the  plan  of  enforcing  boundary  conditions  only  in 
the  net  sense  over  a  portion  of  the  boundary  is  again  neoeseary  in  general. 


¥1.  FORMULATION  IN  TERMS  OF  COMPLEX  VARIABLES 

Formulation  of  the  stress  functions  in  terms  of  analytic  functions  of 
a  oooplex  variable  will  now  be  oarried  out  to  enable  use  of  some  of  the 
elegant  techniques  invented  for  plane  problems.  Since  a  departure  from  the 
Airy  stress  function  is  made,  some  modification  of  the  conventional  oooplex 
variable  representation  will  be  required.  The  two  orders  of  solution, 
n  ■  0  and  n  ■  1,  will  now  be  formulated  in  terms  of  the  oomplex  variable, 

£  p  f  ♦  if). 

The  biharmonic  functions  F,#,(f,ij)  and  G* 1  *  (^,  17)  defined  in  the  pre¬ 
ceding  seotions  can  be  represented  as 

Fm(f,7j)  -  lm|r  jV*’(C)dC  ♦  jV#,U)dt|,  (35) 

-  Ira  If  /Vn(£)d£  ♦  /Vl‘(C)d(;|,  (38) 


where  (£)  and  (£)  are  analytic  functions  of  £,  "Im"  denotes  the 

imaginary  part  of  the  bracketed  expression,  and  barB  denote  the  oomplex 
conjugates.  Equations  11  and  35  yield  the  relations 


(•1 

ae 


(0) 
*  % 


< s>  (si  _  isi  r.  . , isi  -.,isi  ,isn 

"i  -  ~ziTev  ■  -H(3  -  *v)+i  ♦<■*11  +  *1:  J' 

< s  1  r  lei  <»n 

*  %  J* 


(37) 


Similarly,  from  (30)  and  (36)  along  with  (35), 

in  in  r.m  -nil  .  .  »Tm  , (si  -ml 

ai  *  %  ■  J  -  (i/2R)[to{  ♦  -  </»c  -  \  J, 

m  m  ,  rlm  ,  r.m  .ml  ry.m  .i*»l/_ 

,  -  2i T(v  -  -2^(3  -  4K)d>c  ♦  W>cl  t  0t{  J  ♦  ^  J/R 

(C  ♦  t)[(3  -  dv)^*'  ♦  ♦  '/'u'j/SR  *  (1  -  ?’*')/R, 

III  .  Ill  .  in,  ,,  .r.TH>  ,«•» 

o#  ■  v[<r(  +  <rn  )-  (1  ♦  v)[j#c  ♦  W>{  ♦  ♦  ^{  J/2R. 


ill 

-  % 


(38) 
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Curvilinear  coordinates  can  now  be  introduced  by  introducing  appropriate 
oonfaraal  napping  functions  of  a  complex  variable.  Consider  an  auxiliary 
oomplex  plane,  the  y-plane,  where  y  •  e  *  iS.  Let 

C  -  o>(y)  (30) 

denote  a  conformal  mapping  of  a  convenient  region  in  the  y-plane  into  the 
region  under  consideration  in  the  (-plane.  Then  the  analyticity  of  ^,,,(() 
and  \p*M)  (()  ie  preserved  when  considered  as  functions  of  y  due  to  the 
analytioity  of  u(y).  The  argument  is  well-known  in  the  conventional  oomplex 
variable  formulation  and  the  general  details  will  not  be  reproduced  here. 

For  reference  we  shall  list  several  of  the  more  useful  relatione.  If 
on  the  arc  A'B',  we  denote 


d(  •  ds  e 


la 


ds  e 


1 


(40) 


then, 

T1>  +  1  ■  -®  [(3-  4v)<£t  ♦  e 


-*i», 

(<t>,  * 


r.(tl 
♦  L<Pr  r  ♦ 


_m  in  i«r, 

Tv  41  r(3  -  4v)d>{  ♦  e 


-tia  111  -111.  -.Ill 

\<f>r  *  Vr  )  ♦  .  * 


"tt 

111 

"It 


*u  J* 

C3 b 


.  I  • I  ,  I# I .  .  -la,„_<»»  —  It)  I  e I  _l#l.._ 

-  (£  ♦  Z)(tv  *  i  )/4R  -  e  (C?{  ♦  $*><  ♦  )/4R 


)  ♦  w>{ 

nil 

♦  ^  J 

|/2B, 

(42) 

1  tu'*’ 

ml 
**«  j 

i::. 

(43) 

.-**  «4(i  -  v)*'11 

(n/2)  R-t^(3  -  4 v)  fC<t>,$>dC  -  4(1  -  v)i<f>19' -  2(l  -  2v)Z^*} 

,1*1  -...-.1*1  ,1*1, 1»*  ,  , 
♦  ap  -  (C  ♦  D(T<£t  ♦  )|4,,  (44) 

20(u  ♦  i  w  )  ■  (3  -  4v)(d>  -  4>  )  -  •  (45) 
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VZI.  SIllfARY 


A  plan  of  suooesaiTa  approximations  for  handling  the  equations  oT  the 
three-dimensional  elasticity  problem  for  oylindrioal  shells  with  regions 
of  rapidly  varying  thickness  is  proposed  along  the  oonoeptual  lines  of 
0.  Owner.  It  is  shown  that  for  a  restricted  but  useful  class  of  loadings, 
the  utilisation  of  plane  biharmonio  functions  can  be  made;  thus,  many  of 
the  tools  of  complex  variable  theory  oan  be  introduced. 

The  argument  carried  out  in  the  preceding  sections  points  up  the  re¬ 
quirement  for  relaxing  the  precise  loading  conditions  over  portions  of  the 
shell.  This  matter  is  common  to  all  approximate  theories  of  this  type, 
leading  to  the  well-known  separation  of  solution  into  the  so-oalled 
interior  solutions  and  boundary  layer  effeots. 

Finally,  although  the  generally  more  difficult  three-dimensional 
problem  has  been  reduced  to  the  determination  of  plane  biharmonio  funotions, 
the  latter  matter  oan  provide  considerable  difficulties  in  itself. 
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